Conductance through a single impurity in the metallic zigzag carbon nanotube 
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We investigate transport through a single impurity in metallic zigzag carbon nanotube and find 
the conductance sensitively depends on the impurity strength and the bias voltage. It is rather 
interesting that interplay between the current-carrying scattering state and evanescent modes leads 
to rich phenomena including resonant backward scattering, perfect tunneling and charge accumu- 
lations. In addition, we also find a dual relation between the backscattered conductance and the 
charge accumulation. At the end, relevance to the experiments is discussed. 



Single-wall carbon nanotubes (CNTs) [T] are graphenc 
sheet [2] in roll-up geometry with large aspect ratio, mak- 
ing them naturally-made quantum wires with interesting 
transport properties. Since the conducting channels in 
metallic CNT are limited, defects and impurities [3j [4] are 
important for realistic considerations and may play a sig- 
nificant role in transport at nanoscale. One of the recent 
breakthroughs is the technique of covalent attachments [5 
to CNTs, making the fabrication of single-molecule junc- 
tions possible. It is rather remarkable that one can use 
the conductance to monitor and manipulate the cova- 
lent attachments to the targeted CNT. On the theo- 
retical side, vacancy defect or doped boron/nitrogen in 
CNTs [BJ |7] show conductance suppression due to en- 
hanced backward scattering. 

Inspired by these experimental and theoretical stud- 
ies, we investigate the differential conductance through a 
single impurity in the metallic zigzag CNT. Note that 
there are many types of defects including pentagon- 
heptagon pairs, vacancies and chemisorption adatoms[8J. 
For quantitative comparison with experiments, the first- 
principles approaches are more appropriate. On the other 
hand, to reveal the close connections between the conduc- 
tance, the transmission matrix and the current-induced 
charge accumulation, a phcnomcnological approach may 
prove to be useful. As long as the impurity potential 
is short-ranged, the detail profile is not crucially im- 
portant and can be captured by the strength of the 
pseudo-potential, which is taken to be a delta-function 
for simplicity. These simplifications allow us to com- 
pute the scattering states in the metallic CNT and con- 
struct the transmission matrix exactly. Making use of 
the Landauer-Biittiker formula, we can compute the dif- 
ferential conductance at different energies. 

At first glance, this may seem to be a simple calculation 
leading to trivial results - one expects the conductance 
is close to perfect for weak impurity strength, while it 
should get strongly suppressed for a strong impurity po- 
tential. This naive expectation is not correct since impu- 
rities in honeycomb lattice cause strong resonant states 
and non-trival evanescent modes. As a result, the con- 
ductance sensitively depends on the strength of the impu- 



rity potential and the bias voltage. In the absence of the 
impurity, the two conducting channels in metallic zigzag 
CNT deliver the perfect conductance G = 2Gq, where 
Go = 2e 2 //i is the quantum conductance. In some en- 
ergy regime, we found the backward scattering is greatly 
enhanced so that one of the conducting channel is per- 
fectly blocked, leading to resonant backward scattering. 
On the other hand, in some other energy regime, the per- 
fect conductance is recovered, signaling the conspiracy 
from the evanescent modes to make the impurity poten- 
tial transparent. In the following, we would elaborate on 
the detail derivations that provides a firm ground for the 
rich physics of the conductance in CNTs. 

In the following calculations, we concentrate on the 
metallic zigzag CNT that is well described by the tight- 
binding model, 

H=-t^2 (clcj + ctci) + F cjc , (1) 

(id) 

where t is the hopping amplitude between nearest neigh- 
bors and Vo is the strength of the impurity potential 
located at the origin. Since no spin-flip scattering is 
present, the spin just gives rise to a factor of two and 
is suppressed in the above expression. 

The above Hamiltonian can be solved by the general- 
ized Bloch theorem |5| where the momentum is allowed 
to be complex due to the presence of the impurity po- 
tential. In the following, we sketch how this is done 
in metallic zigzag CNT. In the absence of the impurity, 
the transverse momentum is a good quantum number, 
k m = 2irm/N with m = 1,2, ...,N and N is the num- 
ber of unit cells around the CNT. Performing a partial 
Fourier transformation for the wave function in the trans- 
verse direction, <Z>(x,y) = J2k ®k(x)e iky , the CNT is de- 
composed into N decoupled effective ID chains with defi- 
nite transverse momentum k m . In general, the wave func- 
tion after partial Fourier transformation take the form, 
~ 4>kZ x ■ For \z\ = 1, these solutions correspond 
to the familiar plane waves. However, in the presence of 
the impurity potential, \z\ ^ 1 solutions are also allowed 
and need to be taken into consideration. For given en- 
ergy e and momentum k m , z satisfies the characteristic 
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FIG. 1: (Color online) (a) Conductance versus energy for 
(9,0)-zigzag CNT with different impurity strength Vo/t = 
f , 5, 50 shown in blue (dark), red (gray) and green (light gray) 
respectively, (b) The energy-dependent conductance (upper 
curve) and the charge density on the impurity site (lower 
curve) with Vo/t = 5 are plotted together for comparison. 



equation, 



tt, 



(2) 



where t m = It cos(fc m /2). By forming appropriate linear 
combinations of all possible solutions of z, the scatter- 
ing states can be constructed analytically and thus the 
transmission matrix is obtained. Since there are only two 
scattering channels in a metallic zigzag CNT, the 4x4 
5-matrix is 



S(e) 



r(e) f(e) 
t(e) r'(e) 



(3) 



where r, r' , t, t' are 2x2 matrices. Finally, the Landauer- 
Biittiker formula gives the conductance as tracing over 
the square of the transmission matrix, 



G(e) = G Tr [^(e)t(e)] 



(4) 



The transmission matrix can be derived analytically, in- 
volving complicated expressions which can be evaluated 
by numerical means. 

In Fig. [T] the numerical results are shown. First of 
all, since the impurity potential only exists at the origin, 




FIG. 2: (Color online) Charge accumulation at perfect block- 
ing G = Go in the (9,0)-zigzag CNT with Vo/t = 1. The 
radius of the shaded area represents the charge density with 
the defect site labeled in green (lighter gray). 



one can always form an appropriate linear combination of 
the two conducting channels so that one of them remains 
intact from scattering. This ensures that the conduc- 
tance is always large than Go- For weak impurity po- 
tential Vo/t = 1, the conductance is close to the perfect 
value G ~ 2Go at most energies. However, there exist 
two sharp dips down to exactly G = Gq where one of 
the channel is completely blocked. We found the perfect 
blocking is directly tied up with the resonant backward 
scattering with the evanescent modes. Though the con- 
ductance merely comes from the scattering states, the 
evanescent modes are also induced during the impurity 
scattering process. In Fig. [2] the charge accumulation 
at perfect blocking energy is shown with relatively weak 
impurity strength Vo/t — 1. It is remarkable that such a 
large charge accumulation can be created by a weak im- 
purity potential on one site. By increasing the impurity 
strength, the resonant backward scattering is broadened 
but the perfect blocking G = Go remains. Finally, strong 
impurity strength brings the system to the unitary limit 
where the particle-hole symmetry is restored. 

Another interesting feature is the tunneling conduc- 
tance becomes perfect at some energies. This may not 
look very surprising for weak impurity strength. How- 
ever, even for Vo/t = 50, the strong impurity potential is 
still transparent at some particular energies as shown in 
Fig. [T] How can the scattering states ignore the strong 
impurity scattering and maintain the conducting chan- 
nel perfect? To understand this phenomena, we plot the 
conductance and the particle density on the impurity site 
together in Fig. [ijb) for comparison. When the den- 
sity on the impurity site is large, the conductance has 
a dip, while when the density is zero, the conductance 
is perfect. This correlation can be understood from the 
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FIG. 3: Dimensionless backscattered conductance gb = 
Gb I Go and sin 2 r/t versus energy plotted in solid and dashed 
lines. The identical match of these two curves agrees with the 
analytic proof shown in the text. 



Lippmann-Schwingcr equation - with vanishing density 
on the defect site, the scattering state remains the same 
as the unperturbed wave function and thus the impurity 
potential looks transparent. 

The single-impurity scattering consider here is the dual 
junction for the tunneling setup through a disconnected 
quantum dot. In the tunneling junction, the Friedcl sum 
rule relates the conductance to the phase of the transmis- 
sion matrix r\ t and the G = Gosm 2 rj t = Gq sin 2 (7rAiV), 
where r/ t is the phase of det t(e) and AN is the change of 
the particle number. It is rather interesting that a dual 
relation for backscattered conductance is found, 



Gb — 2Gq — G = Go sin rjt, 



(5) 



for the single-impurity scattering as well. We checked 
the relation is valid for our numerical results as shown in 
Fig. [3] According to Friedel sum rule, AN = r] t /n. Thus, 
we come to a useful relation Gb = Gq sin 2 (7rAA r ) that es- 
tablishes the connection between the backscattered con- 
ductance and the charge accumulation. 

To show the relation between Gb and rjt, we start from 
the Dyson equation for the Green's function, 

G wx ,(u>) = G^M + E G° XXi (lo)T XiX2 G° X2X ,(lo), (6) 



where G and G° are the Green's functions with and with- 
out the impurity potential and T xx i is the T-matrix in 
coordinate representation. The non-interacting Green's 
function can be computed by contour integral, G xx , (u>) = 
— mp(<jj)e lk " ][X ~ x ', where k u — e _1 (a; + p) is the mag- 
nitude of the momentum and p(oS) is the density of 
states. Note that the tunneling amplitude t(u>) — 
G XX '(oj)/G xx ,(oj) is the ratio of the Green's functions in 
the asymptotic limit, x — > 00 and x' — > —00, where the 
Dyson equation takes the simple form after some algebra, 
G xx >(uj) = G xx ,(u>)[1 — mp(uS)T(pj)]. By comparison, we 



obtain the key relation between the tunneling amplitude 
and the T— matrix, t(w) = 1 — mp{ui)T{ui). Making use 
of the relation 5=1 — 2-xipT and the unitarity S = e 2lS , 
the tunneling amplitude is t(ui) — cosSe lS . It is clear 
that Gb = Go(l — cos 2 6) and r/t — S so that the relation 
is established. 

Our findings here echo with the recent experiment [10 
that the defect in CNT can turn opaque in the Coulomb 
blockade regime or transparent at other bias voltages. 
Furthermore, there exist controlled methods to fabricate 
the defects 11] systematically, or manipulate the effective 
impurity strength[12 by a voltage pulse from the tip of 
the atomic force microscope on CNT. Therefore, the res- 
onant backward scattering, perfect tunneling and charge 
accumulations discussed here are relevant and can be ob- 
served in experiments on transport through the single 
defect in CNT. 
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